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Abstract 

The simplest version of the Duflo-Zuker mass model (due entirely to the late Jean Duflo) is 
described by following step by step the published computer code. The model contains six macro- 
scopic monopole terms leading asymptotically to a Liquid Drop form, three microscopic terms 
supposed to mock configuration mixing (multipole) corrections to the monopole shell effects, and 
one term in charge of detecting deformed nuclei and calculating their masses. A careful analysis 
of the model suggests a program of future developments that includes a complementary approach 
to masses based on an independently determined monopole Hamiltonian, a better description of 
deformations and specific suggestions for the treatment of three body forces. 
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I. INTRODUCTION 



Masses are a fundamental property of nuclei, whose accurate knowledge is important for 
a large number of processes in nuclear physics, in particular in astrophysical phenomena 
and for a variety of applications in other areas, from elementary particle physics to a precise 
determination of the kilogram. Though much progress has been made in measuring the 
masses of exotic nuclei (see for example Ref.p] and references therein), theoretical models are 
still necessary to predict them in regions far from stability [3] . Advances in the calculation of 
atomic masses have been hampered by the absence of a full theory of the nuclear interaction 
and by the difficulties inherent to quantum many-body calculations. There has been much 
work in developing mass formulas with either microscopic and macroscopic input or within 
a fully microscopic framework. 

Soon after the latest compilation of nuclear masses AME03 was published, a compar- 
ison between the predictions of a large set of mass models, which were fitted to describe 

n 

the nuclear masses included in AME95 [5|, for the more than 300 new masses included in 
AME03 was presented \^. It became evident that the predictions of the Dufio-Zuker (DZ) 
model 6| were outstanding. Recent tests of the predictive power of nuclear models [.7ii| con- 
firm the ability of DZ to make stable predictions which are more accurate than those offered 
by other models. 

The model has many parameters: 28 in the original version, tabulated by Audi and 
used in the comparisons [3J . The code is unpublished but a version with up to 33 parameters 
circulates and it is now available in js]. Another code, with only 10 parameters (DZIO) 8] 
leads to good though less spectacular results but embodies the essence of the model. 

Our task is to examine in full detail what DZIO does and as we go point to possible ways 
to improve on it, based on what has been learned in recent years. To anticipate on what 
follows: Of the ten terms enumerated in the abstract the three — crucial ones — supposed 
to represent configuration mixing are anomalous in the sense that they contain three body 
effects of unknown origin and they scale with the total number of particles A instead of A^^^ 
as shell effects should. They will be examined in special detail. 

In Section [TTl the basic elements of the DZIO model are described. Section UTT] deals with 
the evolution of the results as the ten parameters of the fits are switched on; the origin 
of the anomalous scaling becomes apparent. Predictive power and stability of the model 
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are examined in Section IIVI An analysis of the action of the anomalous terms follows in 
Section |Vl A full explanation of their nature is given in Section I VI I by comparing the 
macroscopic terms with an independently determined monopole Hamiltonian. Section I VII I 
is devoted to the need to introduce three body forces. Section IVIIIl is the conclusion and an 
homage to Jean Dufio. 

II. BASIC ELEMENTS 

To illustrate what a model of nuclear masses is supposed to do consider Fig{T] where 
experimental binding energies of even-even nuclei are subtracted from an improved version of 
the Bethe-Weizsacker liquid drop (LD) form, given in Eq. ([1]) for Z protons, and neutrons, 
with mass number A = N ^ Z and isospin T = \N — Z\ j1. The asymmetry term is 
generalized to include Wigner and surface contributions. Only even-even nuclei are shown 
as they contain all the basic information. The pairing term is omitted leading to (mostly) 
positive definite shell effects. 
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FIG. 1: Shell effects (BE(exp)-E(LD)) along isotope and isotone lines (latter displaced by -14 
MeV). Only even-even nuclei are shown. 

E{LD) = 15.5A - 17.8^2/^ - 28.6 — ^ + 40.2 — j^j^ IV^' 

Both graphs contain exactly the same information. Four remarks (injunctions): 
• LD Any model must contain the LD either explicitly or asymptotically. 



• EI The only observed closures worth considering are the "extruder-intruder" (EI) ones 
at N,Z= (14), 28, 50, 82 and 126, in which the subshell of largest angular momentum 
is "extruded" from harmonic oscillator (HO) shell of principal quantum number p and 
becomes the intruder in the p — 1 HO shell. 

• SPH-DEF Flat patterns at the bottom of Figil], as in the regions 90 < N < 110 
and 60 < Z < 75, correspond to well deformed nuclei. Any model must treat them 
specifically. 

• A-*-/^. Shell effects — defined by subtracting any contribution to the total binding 
energy from its asymptotic form — should scale as A^^^. 

These injunctions were incorporated in steps. First came the purely phenomenological 
work of Duflo j^. It was based on simple algebraic forms that fitted masses with RMSD 



of some 350 keV. However, it extrapolated poorly. In a companion paper [10| it was shown 
that the algebraic simplicity could be reconciled with a microscopic derivation provided the 
microscopic shell effects be separated from the macroscopic (LD) form. Then the observed 
closures were taken to define conventional shell model spaces and it was shown that config- 
uration mixing would lead to simple quadratic and quartic forms in the number of valence 



particles lOj. The trouble was that the observed closures had to be put by hand. The 



breakthrough came by incorporating hints from ll| showing that from realistic interactions 



one could extract a "master term" containing both the leading LD bulk energy and strong 
harmonic oscillator (HO) shell effects. It was left to shift the HO closures to the observed 

n 

ones and then use the shell model forms derived in [10| to mock configuration mixing. In a 
recent paper [12j] the DZ strategy is presented in some detail, stressing in particular that it is 
not a "mass formula" but a functional of the orbital occupancies and explaining its success 
in dealing with well deformed nuclei. Our aim here is to show how to incorporate the last in- 
junction above, by unearthing the reasons why DZ fails to respect it, which we (cryptically) 
anticipate with the benefit of hindsight: DZ assumed that asymptotic behaviour involves 
an LD form which ignores A^^^ terms that turned out to be present and treated as genuine 
shell effects. Next section starts dealing with the problem. 



4 



A. The Master term 



In DZ [6] it was assumed (guessed) that reahstic two body interactions generate two 
collective terms solely responsible for the leading LD contributions, of the form 



Ma 



huj 



Mt 



hjjj 



(2) 



p / 'J \ p V p 

where hw is the HO frequency, hojo is left as a free parameter. Dp = {p + + 2) is the 
degeneracy of the major HO shell of principal quantum number p, nip = Up + Zp, tp = Up — Zp, 
where Up, Zp are number operators for neutrons and protons, respectively. 

To obtain asymptotic estimates for Ma we assume at first N = Z and, following Ref. [3], 
sum up to the closed Fermi shell pf, which we associate to the total number of particles A 
through (use p^^^ = p{jp — 1)(]9 — 2); =^ leads to, implies; ~ approximately; .'. therefore) 



Pf 



A = ^ mp = ^ 2Dp 

p p=0 



2(j9/ + 3)(3) 2{pf + 2f 



{pj + 2)^{UY/' 



Next we estimate hw by relating it to the observed square radius 



13| 



5^mp(p + 3/2) 



AMnucl^^ 



^1/3 

ipf + 2):. /icu = 35.59-^MeV. 



(3) 



(4) 



where Mnuci is the nucleon mass. The 35.59 coefficient follows form an accurate estimate of 
the nuclear radii of magic and semimagic nuclei [l^ , which may serve as a reminder of the 
need to incorporate an isospin dependence. 



- A^-' ( .943 - 0.4-A_ _ o.34(A)2 ) ^(i.04M). 



t = N -Z. 



(5) 



In what follows, whenever convenient, we shall replace hu/huo or hui by the scaling factor 
1/p suggested by Eq where 

2" 



0.5 



(6) 



retaining the DZIO choice for (r^) rather than the more sophisticated Eq. ([5]) 
Approximating ~ p + 3/2, leads to 



Pf 



Pf 



J2^y^p^Pf{pf + A)^iPf + 2f 



(7) 
(8) 
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Similarly, for N and Z separately (= equivalent) 

where pj^, pf^^ are the neutron and proton Fermi levels. As 

Ma = Mz + Mp + 2^/M^^/M^, Mt = Mz + Mp -2>/M^^/M^, (12) 
the leading asymptotic (x) estimates in t/A become 

To go beyond leading order is delicate, but combining Eqs. ([3]) and ([7]) we find that A*^^ 
goes as pj + 8pj + 23pj + . . . against pj + 8pj + IGp'j for rripDp ^^'^)'^, which points to 
a substantial A^^^ pj + 2) contribution to M4 but a vanishing one in A'^^^. In other 
words — as it stands — the master term has no contribution to the LD surface energy. To 
look for its possible origin we examine the microscopic derivation of Ma- 



Z^^J^, (9) 

(PZ^+2r ^ (3^ 
4p ^ 4AV3 ' 



(11) 



B. Origin of the master terms. Scaling 

Let us rely on the general factorization property for arbitrary quadratic forms [ll| 

J2 V.a^. ■Q^ = J2eJY1 ^^■^'^m E ^/^^Z^, J , (14) 

x,a n \ k /3 / 



which we specialize to the monopole part of the Hamiltonian (detailed discussion in |15|). 
in which case V are symmetric matrices, diagonalized by unitary transformations U and Q 
stand for isoscalar (number m), or isovector (isospin T) operators. To fix ideas consider 
the result of diagonalizing the isoscalar monopole interaction for the first 8 major oscillator 



shells for the chiral N3L0 interaction 



16[ smoothed by the V^owfc procedure [l7|. There 



are 36 subshells and as many eigenvalues. One of them turns out to be strongly dominant. 




principal quantum number p 

FIG. 2: The isoscalar master eigenvector (boosted by an arbitrary factor) for the N3L0 interaction. 
According to Eq. (jlSp only the average term in Up contributes to the closed major shells. 



Within a very good approximation its value is proportional to ^ = 1/p and its eigenvector 
is independent of it. The corresponding factor in Eq. (fT^ can be split as 

Y^rrikUk = nipUp + Y,^k{Uk - Up), Up = LkMI^^ Y.^^ = SZ^^Jfc + 1) = Dp, (15) 

k k ^ k k 

where jk is the angular momentum of subshell k. Only the term in Up survives when 
the major shell is full. It defines the collective monopole operator in full analogy with its 
multipole counterparts such as pairing and quadrupole and, again, in analogy with them 
15| we expect Up oc Dp (oc proportional to). The calculated Up are rescaled by a factor 



six so as to make them of order unity and two fits are made to the t = patterns, yielding 
Up = A.2hD~p^^'^ , and a variant U^ = AAYDp^^"^ — O.QDp^, which are seen to be almost 
undistinguishable in Figure [2l However, the asymptotic contributions to the master term 
are quite different Ma x 17MA - 20.97^^/3 and x 17.27A - lO.SlA^/s - 5MA^/^. 
They will referred to as "asymptotics" for short, to distinguish them from "shell effects". 
In the first panel of Fig. [3] it is seen that both approximations produce nearly the same 
shell effects that show as parabolic segments bounded by HO closures at = 8, 20, 40, 70, 
112 and 168. A remarkable result, in view of the very different asymptotics. Each segment 
can be represented by the form mp{mp — Dp){Dpp)^^ , where the denominator ensures at 
midshell (m^ = Dp/2) an amplitude Dp/Ap oc A^^'^. To check that this is the correct scaling, 
in the second panel of Fig. |3] the master shell effects of the first panel are subject to an 
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FIG. 3: Master shell effects produced by Ma and for t = N - Z = 0. See text. 

A^^^^ compression. It is seen that the amplitudes become constant beyond N, Z 70. The 
deviations from perfect (uniform) scaling may have consequences that deserve a comment. 

Scaling is an asymptotic notion: p + 2 and, say, p — 3 both scale as A^^^ but are quite 
different for small p. They differ in a "unit scaling", i. e., a. factor that goes asymptotically 
as unity, in our example {p — 3)/{p + 2). The difference between the two approximations to 
Up in Figure [2] are of the same type. Fig. [3] provides another instance of the possible origin of 
unit scalings, which in DZ work are often represented by surface terms: typically an operator 
F affected by some unit scaling is replaced by F(l + a/p). More generally one may consider 
rational functions i. e., quotients of polynomials of the same rank in p (avoiding vanishing 
denominators, as in Eqs.( |T8|) in next section). 



C. The HO-EI transition 



In the DZ implementations the combination M — )■ Ma + Mt is employed, explicitly 

2 / \ 2' 



M = — 
2p 



nir. 



p 



tr. 



P 



(16) 



while the asymmetry terms are represented by 4T(T + 1)/A, (T=|t|/2) with the shell 
effects in Mt simply ignored. This neglect deserves reexamination and here we only note 
the asymptotically correct form in T(T + 1) rather than is both theoretically sound and 
empirically significant. 

The M term raises one of the outstanding problems in nuclear physics: realistic inter- 
actions fail to produce the observed closures pLS], which contradicts the basic tenet of the 
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discipline that they are due to the spin-orbit force: what the famous / ■ s term does is to 
make the orbit with largest angular momentum the lowest in shell p. This term is 

indeed present in the realistic potentials and it suggests the correct magic numbers but it 
does not produce them [l^. The necessary mechanism has to be invented. Figure H] indicates 



P+ 1 



r{p + 1) 




FIG. 4: Harmonic oscillator and extruder-intruder (EI) shells. 



what has to be done: To change the HO closures (at N,Z=8,20,70 . . . ), into the observed 
extruder-intruder (EI) ones at N,Z=28,50,82 and 126. As made clear by Fig. [T] these are 
the only obvious ones. Therefore, the only relevant operators must separate orbit j{p) of 
degeneracy -Dj(p) = 2{p+ 1) from its partners r{p) of degeneracy -Dr-(p) = p{p+ !)• The only 
one body operators that do it properly are 

(17) 

because they vanish at HO closures and therefore give no asymptotic LD contribution. The 
2{p+l) denominator is arbitrary, u and tt stand for neutron and proton orbits, n and z for 
neutron and proton numbers. Duflo invented the following operator 





'puj^ -2n^/ 








[ 2(p+l) J 




[ 2(p+l) J 



p^ + Ap — h ^ p^ — Ap + h 



E 



5^ - 



'7rp 



p^ + Ap —'b 
^p{p + 2) 



D,ip + 2) 
2 - 4p 5 



p 



TVp 



Dp{p + 2) 



S 



p 



(18) 



which represents a sophisticated example of unit scaling discussed at the end of the previous 
section. It leads to the remarkable result in FiglSl where the HO peaks are practically erased 
to give way to EI ones, which, from left to right, correspond to = 74 (Z = 50), = 82, 



N = 106 {Z = 82), N = 126, = 150 {Z = 126) and = 184. The choice of the 
N — Z = 24 is arbitrary. 
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FIG. 5: The evolution from HO (dots) to EI (squares) shell effects for N — Z = 24 even-even nuclei. 
The asymptotics are roughly represented by a simple A term. Heavier marks refer to species whose 
masses have been measured. 



D. The DZIO equations. Macroscopic terms 



The discussion in Section III Bl suggests that volume and surface effects be treated in a 
single term but in DZ they are kept separate (which turns out to be a problem, as will be 
explained in the comments to Fig. [6]). The terms are: 



M 

M + S and — respectively. 
P 



(19) 



The Coulomb term includes the charge radius Tc and a surface correction (See Ref. for 
more sophisticated treatments) 

2" 



= ZiZ-l) + 0.7Q[ZiZ-ir/^] _ ^ ^,/3 



(20) 



The asymmetry term uses the T(T + 1) form, sounder than the usual (A^ — Z)"^ 

4T(T + 1) 



Vt 



(21) 
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The surface asymmetry has a small unconventional part that can be viewed as the only 
purely phenomenological — representative of isovector shell effects, 

4T(T+1) 4T(r-i) 



Vn 



TS 



The Pairing term includes corrections of order 2T/A: 



N 



Vp 



(22) 



even even (2 — 2T/A) / p 

even odd N > Z (1 - 2T/A)/p 

odd even N > Z 1/p 

even odd N < Z 1/p 

odd even N <Z (1 - 2T/A)/p 

odd odd 2r/(Ap) 

Though purely phenomenolo gica l, these corrections reflect pairing effects of Coulomb and 
nuclear isospin breaking origin 19|. They ensure that the four mass sheets are described 
with equal accuracy. The 1/p scaling in y4~^/^ replaces the now deprecated A"^/^. 

The contribution to the binding energy is thus defined by six terms, of which the last 
four are very close to the usual LD form: 



M 

(H^) = ai{M + S)-a2 a^Vc - Vt + Vts + ^6 Vp. (23) 

P 



E. Microscopic sector. Anomalous spherical terms 

The EI spaces defined by the macroscopic (macro) sector are treated as model spaces in 



which to perform Shell Model calculations. The issue is discussed in Ref. 10|. Here we offer 
a compact summary. 

To estimate the energies of an exact wavefunction |0) we write it as an unperturbed part 
|0) acted upon by fc-body scalar correlation operators Ak (the evolution operator exp (iHt) is 
a scalar). The resulting energy has the form given in Eq.f l24p . where we separate a diagonal 
part given by the monopole Hamiltonian and a correlation one involving the multipole 



Hamiltonian H_ 



M 
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15|. Both Hjji and Hm are effective operators adapted for work in 



the valence spaces. Note that Ai cannot contribute since the space does not allow scalar 
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particle-hole excitations (no two orbits have the same j and parity). 



|0) = (1 + 5^i,)|0) ^E= {0\Hm\0) + {0\HmA2\0) 



(24) 



k 




(25) 



In Eg. (12 5 P we use the notations ml, = m^im^ — 1) and = Dy — to write the form of 
the possible contributions in a valence space v of degeneracy D^. The macro term is solely 
responsible for the closed shells at m„ = = 0. The origin of the quadratic term can be 
both monopole (as in Fig. [3]) and multipole (e. g. the pairing interaction in its simplest 
form). We have no argument to ascribe the cubic terms to configuration mixing. At the time 
the DZ model(s) were formulated the possibility that it could be due to a genuine monopole 
three body force was thought to be unlikely; nowadays it has become a certainty that such 
forces are essential. In section IVTl we shall find clues indicating that the quadratic term owes 
much to configuration mixing and that the cubic one must be due to a genuine three body 
force whose determinant role in the HO-EI transition will be established in Section I VI II We 
are left with the quartic correlation term whose form is dictated by the remark that HmA2 
is a four body operator that must vanish at = rn^ = or 1 because only acts at 
these points. 

In Eq. ( 125|) the denominators are chosen so as to ensure correct A^^^ scalings. Unfortu- 
nately, an original error (explained in the discussion of Fig. |6]) leads to the need of anomalous 
scalings for the "spherical" terms that operate in the EI valence spaces. Using u and vr in- 
deces for neutron and proton quantities, they take the form 



■S3 = - 



1 \ni,n^{ni, - rii,) 



+ 



(26) 



P [ D, 




(27) 



leading to contributions to the binding energy 



{Hs) = 0753-03 h 09^4. 

P 



(28) 



F. 



Deformation term 



That the onset of deformation is due to the promotion of four neutrons and four protons 
from an r shell to the next major HO shell is something that can be read from Nilsson 
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diagrams as pointed out in [lO|. The loss of macroscopic (monopole) energy is upset by 
the gain due to the quadrupole force, simulated by a specific quartic operator (which scales 
correctly). Calling generically n' = n — 4, n' = n + 4, its form and contribution to the energy 
are 



D 



3/2 



3/2 



{Hd) = aio c?4. 



(29) 



The idea was very successfully incorporated in 



the spherical description of rotational nuclei 



12 



;he original DZ fits, and is at the origin of 
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20|, by now firmly established, which 



goes well beyond Eq. (!29|) and points to its limitations. The valence space involves two 
contiguous major HO shells in which good approximations close to Elliott's SU3 symmetry 
operate. Eq. ( |29l) simulates well the mechanism at the beginning of the deformed regions 
by staying in the EI space where the r particles represent the lower shell and the j orbit 
the upper one. It fails as soon as the r orbits are full. Fig. [12] in Section |V] will provide an 
example of the problem. 



III. THE FITS 

Two calculations are made for each nucleus. Both include the macroscopic contribution 
plus either the anomalous spherical terms or the deformed one. 

The coefficients are varied to minimize the root mean square deviation (RMSD) be- 
tween the predicted binding energies BEth{N, Z) and the experimental ones BE^xpiN, Z), 
reported in j4|, modified so as to include more realistically the electron binding energies as 
explained in Appendix A of Lunney, Pearson and Thibault jsj. For each nucleus a spherical 
and a deformed calculation are made, and the one with largest binding is selected: 



BEth = {Hm) + {Hs) if Z < 50 
BEth = {Hm) + max{{H,),{Hd)) if Z > 50 



(30) 



RMSD 



E[5Ee.p(Ar, Z)-BE,^{N, Z)] 



2 ^ 1/2 



nucl 



(31) 



Nnud is the number of nuclei for N, Z > 8. The minimization procedure uses the routine 



Minuit 



|2l|. 
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A. Macroscopic sector 



The heavy task of the macroscopic sector is to ensure asymptotically a LD form and at 
the same time generate shell effects that move the HO closures to EI ones. As the Coulomb, 
asymmetry and pairing terms are already of LD type — within minor provisos — the task of 
generating shell effects falls on the first two terms in Eq. (!23|) in charge of bulk volume and 
surface contributions, as we have seen in Fig. |5l Table [J follows the evolution of the fits as the 
coefficients are turned on. If we judge by these numbers, DZIO does hardly better than any 

TABLE I: The macroscopic coefficients, with their associated mean and RMS errors. 



Op 






M + S 


ai 


9.3980 


5.2535 


4.6971 


16.6714 


17.3653 


17.5337 


M/p 




0.0000 


-22.9448 


-24.6563 


11.8321 


14.7737 


15.4380 


Vc 


03 


0.0000 


0.0000 


-0.0405 


0.6680 


0.6870 


0.6946 


Vt 


04 


0.0000 


0.0000 


0.0000 


26.1232 


35.6940 


36.1628 


Vts 


05 


0.0000 


0.0000 


0.0000 


0.0000 


47.0442 


48.4240 


Vp 




0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


5.0943 


RMSD 


64.062 


29.190 


29.063 


4.232 


2.934 


2.852 


mean 


21.032 


-2.222 


-2.359 


-0.163 


-0.043 


-0.053 



standard LD fit. However, this is not an LD fit and surprises are in store. Fig. |6] displays the 
differences between the experimental and theoretical binding energies, calculated using only 
the macroscopic terms, Eq. ( |23|) . The pattern is both disconcerting and reassuring. It blows 
up at around ^ 80 and there is no trace of the expected A^/'^ scaling. Things go much 
faster and the same time become nicely systematic. The scaling problem (the "original 
error" mentioned in Section IIIE|) can be traced to Section IIIBI where it was shown that 
different variants of the master term have the same shell effects but different asymptotics; 
the leading M + S term in Eq. ( 1T9|) has the correct asymptotics in A, no A^^^ surface term 
and a A^/^ part whose effects are sizeable (about 100 MeV around A = 100). The M/p 
surface contribution in Eq. ( [T9|) — which also contains HO shell effects — is left in charge of 
correcting the A^^^ "drift". An artificial procedure that explains the need of anomalous 
scalings but does not hide the systematic behavior they have to simulate. 
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FIG. 6: Differences between tlie binding energies predicted by DZIO macro, Eq. ()23p and tlie 
experimental ones. Even-even nuclei (RMSD=2.86 MeV). Lines join points at constant t = N — Z 

B. The microscopic sector 

In Table HT] the full set of coefficients Oj, with their associated mean and RMS errors, is 
presented. All fits include the six macroscopic terms, plus some, or all, the microscopic ones, 
ay to aio- It is clear that the three spherical terms can provide a reasonable good fit, with 
an RMSD of 0.72 MeV, but the three terms must be present, acting together. With only 
two of them the RMS error is larger than 2 MeV, as is the case if only the term associated 
with "deformed" nuclei is employed. When the 10 terms are active Eq. fl30|) fits the AME03 
set with an RMSD of 0.55 MeV. 

From the the discussion of Fig. [6] we understand both why anomalous scalings become 
necessary and why the terms of Eq. ( I28l) manage to be so efficient. The factor 2v^+v^ in 
Eq. ( 127|) . though atrocious, is a symptom of the difficulty of accommodating Fig. |6]with a 
single scaling. Eliminating it costs only some 50 keV. On the other hand, forcing the A^^^ 
scaling leads to losses of some 200 keV. We shall examine closely the action of these terms 
in Section IVl 

C. The role of deformation 

From the 2149 nuclei with masses reported in AME03 and N,Z > 8, 1827 are found to 
be spherical and 322 deformed. Figure [7] shows the DZIO deformed and spherical binding 

in 



TABLE II: The full set of coefficients, with their associated mean and RMS errors. 



Op 






M + S 


ai 


17.492 


17.542 


17.769 


17.778 


17.770 


17.766 


M/p 


0-2 


15.284 


15.507 


16.258 


16.355 


16.210 


16.314 


Vc 


«3 


0.693 


0.694 


0.708 


0.708 


0.707 


0.707 


Vt 


04 


35.513 


35.721 


38.354 


37.480 


38.080 


37.515 


Vts 


05 


45.836 


46.653 


56.734 


53.232 


55.394 


53.351 


Vp 


«6 


5.414 


5.275 


5.5361 


6.373 


5.269 


6.199 


S3 


«7 


0.062 


0.448 


0.000 


0.390 


0.000 


0.478 




«8 


0.000 


2.106 


0.000 


1.763 


0.000 


2.183 


Si 


«9 


0.000 


0.000 


0.0215 


0.025 


0.000 


0.022 


di 


aio 


0.000 


0.000 


0.000 


0.000 


37.568 


41.338 


RMSD 


2.443 


2.293 


2.028 


0.717 


2.280 


0.554 


mean 


-0.040 


-0.043 


-0.035 
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energies subtracted from the experimental ones for four chains of isotopes. The crossings 
signal the onset of deformation, which reproduces perfectly the N=90 transition region. 
Note as a pleasant result the correct inclusion of ^^''Nd among deformed nuclei. It should 
be noted though that DZIO underestimates the number of nuclei involved, as anticipated in 
Section IIIFI and illustrated in Section |Vl 

IV. PREDICTIVE POWER AND STABILITY OF DZIO 

The quality of a fit depends not only on its RMSD but also on its error pattern: The 
closer the latter is to a uniform random distribution the better. As we have seen in Fig. El 
systematic deviations are not necessarily a bad thing. Strong isolated discrepancies are more 
serious. Figs. [8] shows that both effects are present in the full fit. To assess their impact 
we propose to follow some recent work [3] where a number of tests were introduced which 
probe the ability of nuclear mass models to extrapolate, all based on nuclear masses taken 
from ref. |^]. The tests are performed for the 2149 nuclei with > 8, Z > 8 or the 1825 
nuclei with A^ > 28, Z > 28. They are: 



Ifi 



-10 

-12 
-14 
-16 



sph 
del 



2 

1 

-2 
-4 
-6 



-10 
-12 
-14 
-16 

-18 



sph 
del 



95 
N 



100 



105 



-10 
-12 
-14 
-16 

-18 



sph 
del 




c d 
FIG. 7: Transition from spherical to deformed for the following chains of isotopes (a) Neodymium 
(eoNd); (b) Samarium (62Sm); (c) Ytterbium (70 Yb) and (d) Uranium (92U). 

• AME95-03: The subset of nuclei with measured masses in the AME95 compilation |5| 

n 

is fitted. This test was used in ref. [3|] to compare predictions of different models. In 
this work the actual masses used in the fit were taken from AME03, only the set of 
nuclei to be fitted is based on AME95. 

• Border region: Nuclei which are furthest removed from stability are excluded from 
the fit and subsequently predicted by extrapolation. 

• Lead region: Nuclei with mass number A < 160, 170, 180, 190 or 200 are fitted and 
the remaining ones, which always include the region around ^"^^Pb, are predicted by 
extrapolation. 

The number of nuclei with predicted masses ranges from 371 in AME95-03 to 810 for 
A < 160. The RMSD of the 16 fits are summarized in Table lllli In the first line the 
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FIG. 8: Differences between the experimental binding energies and those calculated with the full 
fit. (RMSD= 0.55 MeV) 

RMSD of the fits of the whole set of nuclei are presented. The fit is better when the lighter 
nuclei are excluded. Notice that, while the RMSD of the fits is slightly smaller for the 
subsets AME95 and border, as compared with the full set of nuclear masses, when the 
region around ^osp]^ jg excluded from the fit the RMSD are always larger. Also for the three 
subsets A < 160, 170, 180 the RMSD of the predictions is noticeably larger than the RMSD 
of the fits. The situation changes drastically for A < 190, 200, where the RMSD of the 
predictions are slightly smaller than those of the fits. 

Associated with the 16 fits described above are 16 sets of coefficients {cij}. They allow for a 
statistical analysis of their mean value, their root mean square deviation, and the percentage 
of variation, estimated as 100 a/ \mean\. These numbers are reported in Table HVl In the first 
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TABLE III: RMSD of the different test applied to DZIO. 



Test 


N,Z>8 


N,Z> 28 




fit 


prediction 


fit 


prediction 


Full set 


0.5537 




0.4819 




AME95-03 


0.5076 


0.7823 


0.4469 


0.6276 


Border 


0.4988 


0.8212 


0.4540 


0.6527 


A < 160 


0.5787 


0.9049 


0.4823 


1.2330 


A < 170 


0.5847 


0.7740 


0.4966 


1.2889 


A < 180 


0.5708 


0.7711 


0.4859 


1.0481 


A < 190 


0.5855 


0.5611 


0.5098 


0.5010 


A < 200 


0.5844 


0.5014 


0.5054 


0.4728 



TABLE IV: Coefficients of the DZIO mass formula obtained from the best fit of the full data set, 
their average, dispersion and fractional variation from the 16 fits. 



Coefficient 


ai 02 as 04 05 


full set 
mean 
a 

100 a /mean 


17.766 16.313 0.707 37.514 53.344 
17.773 16.332 0.707 37.327 52.277 
0.011 0.027 0.001 0.333 1.759 
0.06 0.17 0.19 0.89 3.37 


Coefficient 


Og O7 Og O9 Oio 


full set 
mean 
a 

100 a /mean 


6.1985 0.4784 2.1831 0.0216 41.3423 
6.2206 0.4853 2.1992 0.0218 40.6985 
0.1044 0.0426 0.2045 0.0006 1.4593 
1.68 8.77 9.30 2.95 3.59 



line the values obtained for the full fit of masses wit N, Z > 8 are included for comparison. 
The six coefiicients employed in the macroscopic terms vary less than 3%, the microscopic 
coefficients are found to vary up to 10%. In any case, it is clear that the DZIO is very stable, 
with its coefficients varying smoothly and very moderately when the set of data employed 
in the fit is changed. 
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V. ANALYSIS OF THE ANOMALOUS EFFECTS 



The anomalous spherical terms are conceptually unacceptable but phenomenologically 
crucial. To gain some insight into their behavior we choose five families of data with t = 
N — Z =8, 16, 24, 32 and 40, and examine how the macroscopic (macro) patterns are driven 
to reasonable agreement with the data (exp) in the full fit (dzlO). Fig. |9]for t = 8 shows 
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FIG. 9: Effect of anomalous terms for N — Z=8 even-even nuclei referred to LD. 

that macro comes close to the exp pattern but is slightly unbound. The anomalous terms 
bring in the necessary correction. 
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FIG. 10: Effect of anomalous terms for N — Z=16 even-even nuclei referred to LD. 



For t = 16 in Fig. the situation changes dramatically: macro is way off except at the 
closures. Contrary to the t = 8 family where the naive shell model seems valid, most of 
the t = 16 nuclei exhibit vibrational features. In principle they would demand a separate 
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treatment as done for the well deformed species. As such a treatment remains to be found, 
DZ relies on the anomalous terms to restore reasonable agreement. 




FIG. 11: Effect of anomalous terms for N — Z=24 even-even nuclei referred to LD. 

For N — Z = 24 in Fig. [11] the situation is more complicated: macro gives a good 
description near closures at = 82 and = 106, Z = 82 {^^^Fh) but starts deviating 
strongly at A^ = 90 where deformation sets in. The inclusion of in Eq. fl29l) (macro+def 
in the figure) restores agreement with data for A^ = 90, 92, 94 that rapidly deteriorates at 
A^ = 96. The inclusion of the anomalous terms brings back reasonable agreement. This 
provides a clue that supplements what was found for t = 16. 

DZ does an excellent iob at describing the onset of deformation at A^ = 90 (^^^Dy) 

nnn 

through 4n-4p jumps, an idea vindicated by later work [12|, [151, |20|. However, deformed 
regions not only start somewhere, they also end somewhere, in our case at the ^^*Pb weak 
closure. In between, low lying 7 bands indicate the need to go beyond 4n-4p jumps. DZ does 
not include this option and leaves the anomalous "spherical" terms in charge of deformation 
after fails to do it. Note that contrary to the t = 16 "vibrators" the necessary tools to 
go beyond 4n-4p jumps are available, as explained in Section [11 F[ 

Something similar occurs for A^ — Z = 32 in Fig. [12) macro is good at first, then defor- 
mation sets in and is well described by def for a while but then the anomalous terms take 
over in a region that remains well deformed. They even correct nicely what is missed by 
macro at the A^ = 114, Z = 82 closure. The situation becomes baffling for t = 40 in Fig. [13] 
because of the uncanny capacity of the anomalous terms to restore agreement with data 
that appeared compromised by the distorted macroscopic pattern. 

The anomalous terms can deal with practically everything, including deformation. This 
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FIG. 12: Effect of anomalous terms for N — Z=32 even-even nuclei referred to LD. 
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FIG. 13: Effect of anomalous terms for N — Z=40 even-even nuclei referred to LD. 

is achieved by borrowing the sound idea of simple polynomial forms and subverting it by 
introducing unacceptable scalings. Among the many things they do, the anomalous terms 
correct shortcomings of the macroscopic contribution, as seen in the Z = 82 closure for 
t = 32 and in the t = 40 case. We have already traced the scaling flaws of the macroscopic 
terms in the discussion of Fig. [6] but its worth pushing their analysis in a totally different 
direction. 



VI. ANALYSIS OF THE MACROSCOPIC TERM. GEMO 

Rather than attempting to correct directly the scaling flaws of DZIO, we shall compare the 
model with a complementary approach in which it is very simple to force correct asymptotics, 
and has the enormous advantage of being practically parameter free. 



99. 



Ideally, the macroscopic-microscopic se par ation should mirror the monopole-multipole 



separation of the realistic interactions 
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15| . We recall that the monopole part pro- 



vides the natural definition of the unperturbed Hamiltonian as it contains all the number 
and isospin operators and all that is needed for a spherical Hartree Fock variation. Our 
purpose is to examine to what extent the macroscopic DZIO terms are consistent with a 
monopole Hamiltonian derived without any reference to masses; defined in 22| (DZII or 
GEMO (GEneral MOnopole) after the name of the code jsj, under item containing 22 1]) by 
fitting (quite well) all single particle and single hole states on doubly magic nuclei (the cs ± 1 
set). Some details that supplement the original paper can be found in [12|, ll5[. GEMO is 
fully equivalent to a mean field calculation of monopole shell effects. It decouples approx- 
imately the macroscopic part by replacing the master term by the combination Ma — 4:K 
which has no A and A^^^ contributions and produces shell effects similar to those of Fig. [31 
Ma is from Eq. ([2]) and K is the kinetic energy affected by an huo ~ 11 denominator chosen 
to reproduce the correct coefficient for the Ma — 4:K term. 



hjjj 



{Pf + 3f\pf + 2) 



(32) 



Asymptotically we obtain from a fit (Ma - ^K) x -28.04^^/3 - Q.lht^/A + 9.43; and from 
now on GEMO stands for the properly subtracted value. To compare with experiment we 
subject GEMO to a 2.5 contraction (the only free parameter, to be examined later) and we 
add to it E{LD) from Eq.([T]). In the calculation the filling order of the orbits has been 




FIG. 14: Comparison between GEMO and experimental masses for even-even nuclei. RMSD=1.69 
MeV. Lines join points at constant t = N — Z . See text 



kept fixed throughout, according to an / ■ s law with an / ■ / modulation. Much better 
results are expected if the proper procedure is followed by determining variationally the 
filling order for each nucleus. In spite of this shortcoming the results shown in Fig. [H] have 
substantially smaller RMSD than those for DZIO in Fig. IH] (1.69 vs 2.86 MeV. Beware of 
y-scale differences between the figures). However, for ^ 80 both models have very much 
the same RMSD^ 1.2MeV. The difference is entirely due to a flatter GEMO pattern {i. e., 
more symmetric around the origin), consistent with A^/^ scaling. In the lighter region it 
is impossible to identify common systematics, while the heavier one is dominated by cubic 
parabolic trends, e. g in the EI space between = 82 and 126. 

If we now compare the models for the cases studied in last section we find the following; 
at if: = 8 GEMO is better, at t = 16 DZlOmacro is better, at t = 24 and 32 both are 
comparable provided the filling orders are corrected, at t = 40 GEMO is much better. In 
general, the detailed subshell structure of GEMO is a source of unphysical, jumpy, behaviour 
(common to mean field formulations), while DZlOmacro suffers from distortions due to bad 
asymptotics; but both models are telling the same thing about anomalous terms: 

• Bad scaling is due to bad asymptotics. 

• Cubic terms i.e., three body forces are a fact, not a DZIO artifact. In the next section 
we shall offer a proof of the need of genuine three body forces to explain the HO-EI 
transition. 

Let us conclude with a word on the 2.5 contraction. GEMO was fitted to single particle 
and single hole properties on closed shell nuclei: they demand huge shell effects whose 
relative positions are little affected by total binding energies i. e., by correlation effects. 
Once these are switched on they will produce in first approximation parabolas of negative 
slope m(D — m) opposite to those in, say. Fig. [31 of form m{m — D). The contraction 
factor simulates this mechanism. The argument is purely heuristic — we have not specified 
the nature of the correlations — but conceptually valuable as it provides an estimate of their 
(large) magnitude and explains the impossibility for an effective mean field description to 
describe simultaneously masses and single particle or single hole energies on closed shells. 
In DZIO the correlation effect is directly incorporated in a single master term that produces 
both shell effects and the bulk energy of nuclear matter. A non trivial achievement. 
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VII. REALISTIC FORCES AND THE THREE BODY ISSUE 



The DZ strategy amounts to assume perfect potentials and proceed as if we were solving 
the Schrodinger equation. A perfect potential should saturate i.e., produce the correct 
energies at the right radii. It should also produce the EI closures. It does neither. Therefore 
DZ postulates that the necessary mechanisms can be incorporated into some effective master 
terms M4 and and adopts the standard shell model approach of "freezing" hw at the 



observed radii. The present consensus [15|, |18| demands the introduction of three body 



forces, and in recent years considerable effort has been devoted to derive them from chiral 



perturbation theory, with some success as witnessed by a recent paper [23| where abundant 
references will be found. However, much remains to be done at the fundamental level and 
here we shall follow the policy of consulting the data searching for hints about the monopole 
Hamiltonian. Our first task is to identify the form of the necessary three body terms. 

NOTATIONS: in what follows we use X^t = X,(Xt - + X^^) = X{X - 1). 

We start with an isoscalar monopole two body interaction written as Hm = Yls<t ^siVst 
where the sums extend over all orbits st. From now on we omit the Vst matrix elements 
and consider only the m^t operators. Next we examine the action of Hm in a major HO 
shell, and call c the occupied ("core") orbits below the major shell. It will give three types 
of contribution: m^c' = -Dec', ^cf^s = Dciris and rrist, where we have replaced the core 
occupancies by the degeneracies Dc of the corresponding orbits. We are left with a pure 
core contribution that we assume will be taken care of by the master term, an effective single 
particle term in and the genuine two body rrist part. 

The first key step is to adopt the "invariant representation" 2J, and references therein. 



An early — relatively successful but inconclusive-attempt to introduce many body forces] 
to separate the total number operator rup in HO shell p from others that are taken to be 
"orthogonal" to it in the sense that they vanish for nip = and rhp = Dp — nip = i.e., 
for zero particles and zero holes. This is achieved as follows (the orbit 1 can be chosen 
arbitrarily) , 



nis = nip + r^i^ (33) 



rrist ^ ^mf + (m, - l)ril) + li?, (34) 



where 

(1) _ fnis _mi\ D,Di _ _ frhs _mi\ DsDi _ _=(i) , , 

''~\d, dJd, + d,- \Ds dJd, + d,-'^' ^^^^ 

= isF ^ 5P ' ^ j (J.H-A)'^) = = 

In reading these expressions it is convenient to fix ideas through an example. Consider the 
sd shell. Eq. (133|) gives the spectrum of ^''O. The m part belongs to the master term. The 
single particle energies are referred to one of them through T^^^ operators in Eq. f l35|) which 
change sign when particles are turned into holes in '^^Ca. The normalization ensures unit 
splitting between single particle states s and 1. The m'p ' in Eq. (jM]) goes with the master 
term. Then we have a modulation of the single particle energies (rrip — l)r^^'' that makes it 
possible for the splittings in ■^^Ca to be different from those in ^''O. Finally the T^g^ operators 
which according to Eq. fl36|) are particle- hole sjTumetric and vanish for m, m = 0, 1. They 
reproduce the strictly two body contributions to the centroids (average energies) of two 
particle and two hole configurations. The normalizations are such as to produce splittings 
of order one between the centroids. 

The second key step is as metaphysical as physical: to select very few possible operators. 
If too many are truly needed, our approach is doomed. Much of the success of DZ is due 
to the choice of a single Fp^^ = '^f(^p)r{py '^^^^ masses, but a general needs 

some extra freedom. As any one body operator referred to its centroid becomes "strict one 
body" in the sense that it vanishes at both m^, ffip = two classical choices come to mind 
/ ■ s and /■/ = /(/ + 1) — p{p + 3)/2 (they may be written as combinations of F^^'* 22]). For 



the strict two body case the only obvious choice is Fp^'' = '^f(^p-jr{p)- -^^^ simplicity, in what 
follows we refer to a single generic Fp^'^\ 

So far we have considered one major HO shell. In general, particles may move in different 
shells. The possible cross shell operators will be of the form 

■p(i) -pfii-pCi) 

mpTrip,, mpVy, Tyvy . 

Upon introducing three body interactions we shall encounter terms of the type 

rricc'c" = Dcc'c", TiT'cc'Tn's = Dcc'TUg and Dcirist 

which will modify the two body equations ( 133|M|) . The genuine mgtu parts will contain rrip, 
and Fp^'^-* operators plus eventually cubics Fp^"*. By now we can propose a list of possible 
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operators: 



cross 



intra shell : (ai + birUp + CimJ)rW, (as + b2mp)T^^\ a^T^^'^ (37) 
shell : (amp + AmJ + ^mpmp,)Ty , mp(a2rj,'^ + /^ar^^^rj,))), asr^F^?^ (38) 



We should not forget the pure m contributions that we have decided to ascribe to (the three 



body part of) the master term, see [15|, p 437] for some speculation on the subject. 

Everything we have said applies either to purely isoscalar operators or to a neutron proton 
representation, in which case each major shell is split in two with replaced by and Zg- 
Physically it is preferable to work in an isospin representation involving and Tg operators. 
The reason can be understood by referring to Fig. [1] which shows that shell effects are very 
much the same for both fluids (at constant or Z). If we plot along lines of constant A 
or T, the former become the fairly smooth "/3 decay parabolas" while the latter (of which 
we have seen many examples) exhibit most of the shell structure, indicating its basically 
isoscalar character. Let us examine how to proceed. 



A. The origin of the HO-EI transition 

Once we accept that three body effects are necessary, we must find ways to identify 
the relevant operators. As there are several candidates, to select them by fitting masses 
may give non unique answers. The best way to proceed — as far as the HO-EI transition 
is concerned — is to look into spectroscopic information: an apparently impossible task that 
turns out to be very simple. 

The unequivocal sign that a strict two body treatment is not sufficient came from exact 



calculations leading to a J'^ = 1^ ground state of instead of the observed 3"*" [25|, l26| . 
a long standing puzzle in conventional shell model work, where exactly the same problem 
exists in ^^Na. Then it was shown that by changing the Vjj, Vj'j^ centroids of a realistic 
interaction R according to 

v;',{R) ^ v^A^) - i-r K ^^^^ 

one could correct much of the spectroscopic trouble in the p, sd and pf shells. To within 
some fine tunings — and a major difference — this is a time- honored two body prescription {28 1 
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that has become a common feature of effective interactions in the pf shell [15|, sections V, 
VB]. The major difference is that k is a linear function of rrip, which makes the prescription 
three body. But it (still) has a flaw: as it was originally devised to ensure EI closures at '^^Ca 
and ^^Ni where the centroids play no role, they were left out. To include them we note 
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FIG. 15: Spectra of Na and Si. The uncorrected and corrected Bonn C interaction (BC) 29 1 



is compared with USD 



and experiment. 



that the six centroids involved, once cast as three isoscalar and three isovector terms, reduce 
to four by removing the overall m(m — 1) and T(T+ 1) dependences. According to our policy 
of omitting at first the isovector pieces, we are left with the terms in Eq. fl37|) . We ignore 



the last and redo the calculations in 27j. We know from 2^ that (ai + birrip + cimp)r 



plays a role, but it does not solve the problems. Therefore we simulate its influence by 
pushing up the ^3/2 orbit two MeV above its observed value (arbitrary move but of little 



consequence). We are left with (02 + &2'^p)r 



(2) 

j{p)r{p) 



K(m)T 



(2) 

j{p)r{p)' 



which is seen in Fig. 



to work very well: the dismal spectra produced by the original interaction (Bonn C k = 0) 
become adequate in ^^Na (Bonn C k(6) = 0.9) and good, even when compared with USD, 
in 27Si (Bonn C ^(11) = 0.6). There is no way to find a compromise, constant, k. 

A refined approach will demand the inclusion of some other operators, but there are 
reasons to believe that K{m)T^^ will play a (the) fundamental role in the HO-EI transition. 



VIII. HOMAGE TO JEAN DUFLO AND CONCLUSIONS 



In all justice DZIO should be called DIO because it was entirely the work of the late Jean 
Dufio. It was designed to shed some light on the problems of DZ28: 
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a) Too many parameters; b) Obscure origin of EI closures; c) Many surface terms that 
lead to a change of sign at around A — 120; d) Anomalous A scalings. 

At the price of a substantial but acceptable RMSD increase, Duflo solved a, clarified b, 
reduced c to a single three body term and was left with d. He had a genius for regrouping 
and eliminating theoretically plausible contributions and inventing phenomenological ones 
demanded by data. Jean was a magical data manager. He was also aware that DZIO needed 
improvements, and to the last day he worked on them. He did not find a satisfactory 
solution, but from his notes one could guess that he did not worry about anomalous terms 
and concentrated on basic shell formation. Our present work points to the same direction. 
We have explained how the scaling anomaly can be eliminated and demonstrated the need 
of genuine three body forces. It remains to be seen whether these advances can lead to a 
formulation as compact and successful as DIO. 
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